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ABSTRACT 



Aims. In this second paper of the sequel of two papers, we further investigate the problem of molecular cloud (MC) collision. 
Anathpindika (2009) (hereafter paper I) considered highly supersonic cloud collisions and examined the effect of bending and shearing 
instabilities on the shocked gas slab. We now consider moderately supersonic cloud collisions (precollision cloud velocities of order 
1.2 km s" 1 to 2.4 km s" 1 ). 

Methods. In the current paper, we present five SPH simulations of fast head-on and/or off-centre cloud collisions to study the evolution 
of ram pressure confined gas slabs. The relevant thermodynamics in the problem is simplified by adopting a simple barytropic equation 
of state. We explore the parameter space by varying the pre-collision velocity and the temperature of the post collision gas slab. 
Results. The temperature in a pressure compressed gas slab is crucial to its dynamical evolution. The pressure confined gas slabs 
become Jeans unstable if the average sound crossing time, t cr , of putative clumps condensing out of them, is much larger than their 
free fall time, tjj. Self gravitating clumps may spawn multiple/larger N-body star clusters. Warmer gas slabs are unlikely to fragment 
and may end up as diffuse gas clouds. 

Key words. Molecular clouds (MCs) - pressure compressed gas slab - gravitational instabilities - filaments - star formation - SPH 
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1. Introduction 

Giant Molecular Clouds (GMCs) have turbulent and clumpy in- 
teriors. Potential star forming clouds move randomly and have 
a velocity dispersion of the order of a few km s _1 on the small- 
est scales (e.g. Larson 1981; Elmegreen 1997; Elmegreen 2000). 
One extreme of this problem was investigated in paper I, where 
we simulated highly supersonic, both, head-on and off-centre 
cloud collisions. We found that the dynamical behaviour of a 
shock compressed gas slab was dominated by internal shear. 
Hydrodynamical instabilities suppressed the gravitational insta- 
bility and turbulent mixing between slab layers made the slab 
susceptible to the non-linear thin shell instability (NTSI) and the 
Kelvin-Helmholtz (KH) instability. Further, internal shear led to 
dissipation of internal energy via shocking between fluid lay- 
ers. The shocked slab lost thermal support against self gravity 
and collapsed to form an elongated gas body, along the collision 
axis, quite similar to the integral filament in the Orion A region. 

Star clusters and their birth : Certain features of star forma- 
tion are now clear. It has been observationally well established 
that most star formation occurs in clusters (e.g. Pudritz 2002). 
There are two paradigms of star formation, one that of quiescent 
star formation - the cloud - core picture (Hoyle 1953; Whitworth 
et al 1996) and second, the dynamical picture of star formation 
also called, triggered star formation. 

The former picture is elucidated by self gravitating MCs. 
Potential star forming clouds may become self gravitating due 
to internal perturbations or may be compressed either by shock 
waves or galactic density waves. Such theories are encouraged 
by spatial correlation between between star formation and the 
spiral arms of disk galaxies. Gravitational perturbations grow 
in regions which become Toomre unstable. Further, the surface 



density in star forming regions, E gas , is supposedly related to the 
star formation rate (SFR), "Esfr, through the Schmidt-Kennicutt 
power-law (Schmidt 1959; Kennicutt 1989; Kennicutt, 1998). 

However, this paradigm suffers on two counts. First, it in- 
evitably leads to a correlation between the SFR and the den- 
sity wave amplitude. However, there is no corroborative obser- 
vational evidence (e.g. Elmegreen & Elmegreen 1986; Kennicutt 
1989). Second, this paradigm fails to reconcile the observations 
of star formation in those spirals, which have small density wave 
amplitude (e.g. Block et al 1994). This suggests that the density 
waves and disk instabilities are not quintessential to star forma- 
tion. We may therefore have to explore certain aspects of the 
second paradigm, stated above. Tan (2000) for instance, invoked 
cloud collision in the SFR calculations for a galactic disk and 
arrived at a power-law similar to the Schmidt-Kennicutt correla- 
tion. 

Triggered star formation is a self regulatory process. Rich, 
young star clusters have been observationally well studied. The 
orbital dynamics of clusters have enabled observers to estimate 
various physical properties of multiple systems like orbital ec- 
centricities and orbital stability, stellar masses, among other 
properties. Star clusters have reasonably homogeneous environs 
and therefore, are useful regions to study stellar evolution. For 
a more detailed review of the matter refer (Clarke, Bonnell & 
Hillenbrand 2000). Star clusters however, show a wide variation 
in stellar density. For instance, ~96 percent of the IR sources 
in the Orion B region are distributed over four identified star 
forming regions. On the other hand, in the Orion A region, star 
forming cores are distributed along the north-south direction - 
the Integral shaped filament e.g. Lada et al (1991); Meyer & 
Lada (1999) and Mookerjea et al (2000). 
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Birth of star clusters - Numerical studies : Triggered star 
formation is elucidated by dynamical interaction between fluid 
flows within GMCs. Colliding flows produce gas slabs that 
are dynamically unstable (e.g. Klessen, Burkert & Bate 1998; 
Klessen & Burkert 2000). Dense gas shells swept up by expand- 
ing ionising fronts from massive star clusters are also dynam- 
ically unstable and may fragment to produce a number of po- 
tential star forming clouds and filaments (e.g. Wada, Spaans & 
Kim 2000; Dale, Bonnell & Whitworth 2007; Furuya, Kitamura 
& Shinnaga 2008; Lefloch, Cernicharo & Pardo 2008). 

The other model suggested and investigated further in the 
present work, is that of gravitational fragmentation of an external 
pressure confined gas slab. See for e.g. Elmegreen & Elmegreen 
(1978), Whitworth et al (1994), Clarke (1999), and Boyd & 
Whitworth (2005). Such slabs result from clouds colliding with 
relatively small precollision Mach numbers (e.g. Chapman et al 
1992; Bhattal ef a/ 1998). 

Below, we present five simulations, in that, two each are of 
fast head-on and off-centre cloud collisions, respectively. The re- 
maining one is a slow head-on cloud collision which post colli- 
sion, forms to a continuous shock as against a jump shock in the 
remaining four cases. We shall attempt to contrast the findings 
recorded in paper I with those, presented here. The plan of the 
paper is as follows. In §2 below we briefly introduce the numer- 
ical scheme employed for our work and the code used. In §3 we 
describe the initial conditions and list the simulations performed. 
The results are discussed in §4 and we conclude in §5. The col- 
umn density is measured in M pc~ 2 and all the column density 
plots presented here, are plotted on a logarithmic scale since they 
span about six orders of magnitude. The spatial coordinates are 
marked in parsecs. 



2. The Numerical Method & Code 

We use the Smoothed Particle Hydrodynamics (SPH) numerical 
method for our simulations. This scheme treats the fluid under 
investigation as an ensemble of particles. The simulations pre- 
sented here are performed using an extensively tested SPH code, 
DRAGON (Goodwin, Whitworth & Ward-Thompson 2004). It 
uses the Barnes-Hut Octal spatial tree (Barnes & Hut 1986) to 
search nearest neighbours and, evaluate the net force on a SPH 
particle. The critical cell opening angle used by us is, 9 cr j, ~ 0.45 
and each particle has 50+5 neighbours. In order to improve the 
accuracy of force calculations quadrupole moments of remote 
cells are also included. The code employs the multiple particle 
time stepping scheme (Makino 1991). 

An extremely dense agglomeration of particles in the simu- 
lation is replaced by a sink particle. The sink is characterised by 
two physical variables viz. the sink radius, R S mk, and the sink 
density, p S mk, both of which are predefined (Bate, Bonnell & 
Price 1995). We set p,,-^ = 10~ 12 g crrr 3 in all the simulations 
presented here. The sink radius is so chosen that the initial mass 
of a sink particle is comparable to the minimum resolvable mass, 
M min , in the simulation. M min is essentially the same as the Bate- 
Burkert mass (Bate & Burkert 1997). We had adopted an identi- 
cal prescription in paper I (refer §2 therein). In the present work, 
protostellar cores are modelled as sinks. 



3. Initial Conditions & Numerical Experiments 

We model MCs as unconfined Bonnor-Ebert spheres. Thus, there 
is no external pressure on the cloud edges. Individual clouds are 
assembled in the same way as described in paper I. We compro- 



mise with the intercloud medium in these simulations to permit 
larger number of gas particles in each cloud. 

In the first four models tested here, we adopt the same equa- 
tion of state (EOS) as that in paper I albeit, with appropriate 
density switches. Using the temperature jump condition, we cal- 
culate the post shock temperature, T ps , for each model. Behind 
the shock, the gas is allowed to cool. In the first two simulations, 
the gas downstream of the shock is maintained at original precol- 
lision cloud temperature, while in the next two, the gas is cooled 
down to 10 K. The EOS employed is, 
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Here ks and m are respectively, the Boltzmann constant and the 
average mass of a hydrogen atom, T ps is the post shock gas tem- 
perature. In cases 1 and 2, we set Tq = T c u, while in 3 and 4, Tq 
= 10 K. 

Finally, as a consequence of the very low precollision veloc- 
ity (~ 0.4 km s _1 ) of individual clouds in model 5, the post- 
collision shock is weak. Weak shocks in astrophysical envi- 
rons are nearly adiabatic, and so we model the one here, with 
a slightly stiffer EOS given by 



n r-i 



(2) 



where all the symbols have their usual meanings and the adia- 
batic gas constant, y = | for either EOSs. Models tested here 
have been listed in Table 1 above along with their relevant phys- 
ical details. 



4. Discussion 

The question of gravitational instability in cold, semi-infinite 
planar gas slabs has been extensively studied in the past. 
However, most of these analyses were restricted to slabs with 
no external pressure on them. See for instance, McKee (1999) 
and Larson (1985) . On the other hand, stability analysis of 
cold, external pressure confined isothermal slabs by Elmegreen 
& Elmegreen (1978) and Elmegreen (1989) suggests that the 
minimum mass of a fragment condensing out of such a slab, 
increases with the age of that slab as it continues to accrete mat- 
ter. The monotonic growth time of the fastest growing unsta- 
ble mode in such a slab is ~Q.25(Gpi ayer )~^ (Elmegreen 1989), 
where pi ayer is the average density of the post-shock gas slab. 

Expressions for the length of the fastest growing mode, the 
timescale of growth of this mode and the minimum mass of 
the fragment condensing out of the slab have been derived by 
(Whitworth et al 1994) and Boyd & Whitworth (2005). Below 
we discuss the models tested in this work and compare the re- 
sults of models 3 and 4 with the analytic findings of the latter 
authors. 

Fast cloud collision (Tq = T M ; models 1 & 2) 
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Table 1. List of simulations performed with relevant physical details. 
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Serial 


Experimental 






Precollision 


Number of 


Head-on 


To 


No. 


details 






Mach Number(AI) 


particles 






1 


M cU \ = M cld2 


= 50 M Q 




3 


N gas = 120000 


Yes 


54 K 




Rem = Rcui = 


0.8 pc, T M = 


54 K 








2 


M M = M cU2 


= 50 M 




3 


N gas = 120000 


No 


54 K 




Rcld\ = RcM2 = 


0.8 pc, T cU = 


54 K 








3 


M M = M M2 


= 50M o 




3 


N gas = 120000 


Yes 


10 K 




Rem = Rem = 


0.8 pc, T cli = 


54 K 








4 


M c m = M cM2 


= 2000 M Q 




3 


N gas = 120000 


No 


10 K 




Rem = Rcld2 = 


4.5 pc, T M = 


377 K 








5 


MM = M cU2 


= 50M o 




1 


N gas = 120000 


Yes 


54 K 




Rcldl = R-CM2 = 


0.8 pc, T cU = 


54 K 









Note : The numbers on the left hand side below indicate the respective column numbers of the table. 

(2) Physical details of individual clouds in a simulation (Mass, radius and temperature of individual clouds). 

(3) Precollision Mach number (M). 

(4) Number of gas particles (N gas ) used in the simulations. 

(5) Specifies whether the clouds collide head-on or at a finite impact parameter. 

(6) Temperature at which the post shock gas slab is maintained. 
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Fig. 1. A column density plot of the time (measured in Myr) sequence of cloud collision in model 1. Cloud collision results in the 
formation of a pressure confined gas slab, which then re-expands and terminates in a diffuse gas cloud, as can be seen in the last 



Sna ?n mdtieT r^Dtc^^o? lhass 50 M G each, collide head-on with 
each other at a pre-collision velocity of ~1.3 km s~' (vV( = 3). 
Individual clouds are thermally supported against self gravity. 
The post-collision thermodynamics in this model is governed 
by the EOS defined by equation (1) above. The clouds collide 
and form a planar gas slab confined by ram pressure. The time 
sequence of colliding clouds and their evolution after collision 
is shown in Fig.l, which is a column density plot showing the 
clouds and the post-collision slab. We treat the gas slab as ap- 
proximately isothermal and hold it at the original pre-collision 
cloud temperature. Our assumption of the post collision temper- 
ature should not have any adverse bearing on the final outcome 
of the simulation, since the post-shock temperature, T ps , calcu- 



lated using the temperature Jump condition will only be higher, 
making the slab even warmer. 

The post-shock gas slab in this model does not undergo grav- 
itational fragmentation. After formation, it simply re-expands 
and ends up as a diffuse gas cloud ( t — 3.49 Myrs in Fig. 1). 
Unlike the gas slabs reported in paper I, the one in this model 
does not show any evidence for the bending instability. However, 
gas elements within the slab are non-static and dissipate mechan- 
ical energy via shocking, facilitated by the artificial viscosity. 
The gas slab in this case generally evolves via an interplay be- 
tween gravity and thermal pressure. There are two phases in this 
evolution. In the first, thermal pressure builds up in the slab af- 
ter shocking and then, in the next phase, as we permit the gas to 
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cool, thermal pressure within the slab rapidly diminishes, and the 
slab re-expands. Even then the slab is too warm to become Jeans 
unstable. This seems similar to the re-expansion of the shocked 
slabs reported in paper I. However, unlike those slabs, the one 
in this case does not collapse to form an elongated object along 
the collision axis suggesting that, gravity and thermal pressure 
attain a dynamic equilibrium. 

The gas slab in this case is not sufficiently massive, or in 
other words, it is not sufficiently cold (for the mass that it has) 
so as to support the Jeans instability. According to Whitworth et 
al (1994), the length of the fastest growing unstable mode is 



2a lyer 



-'fastest 



(3) 



'layer 



which is of order a few parsecs for the present case, at least an 
order of magnitude larger than the slab thickness. 

In model 2, we maintain physical conditions identical to 
those in the previous case except, that the clouds now collide, 
off-centre. The impact parameter, b , is a fourth of the cloud ra- 
dius, Rcu- The time sequence of colliding clouds in this model 
is shown in Fig. 2. The cloud collision results in an oblique 
shocked slab unlike a planar slab in model 1. After formation, 
the slab tumbles about the z-axis and eventually re-expands, as 



in model 1 . This can be seen in the snapshot corresponding to t 
= 1.89MyrsinFig. 2. 

As in model 1, the gas slab in this case also evolves through 
an interplay between gravity and thermal pressure. Contrary 
to intuition, the slabs in either cases evolve on a similar time 
timescale. This suggests that the re-expansion of the oblique slab 
also commences at a similar epoch as the planar slab. In fact, the 
oblique slab expands slightly faster than the planar slab of model 
1 . This is presumably due to axial transfer of angular momentum 
as the outer ends of the oblique slab move outward (see snapshot 
corresponding to t — 1.89 Myrs in Fig. 2). The terminal state in 
this case is similar to that in case 1. 

Equation (3) above, is valid for off-centre cloud collision as 
we have chosen a rather small impact parameter. Lf astest , calcu- 
lated using this equation in the present case is much larger than 
the slab thickness. Hence as in model 1, the gas slab in this case 
as well, does not become Jeans unstable. Such a scenario (as in 
models 1 and 2) will thus, simply lead to formation of diffuse 
clouds in the ISM. These clouds may spawn star formation only 
when they have been sufficiently squeezed. 

Fast cloud collision (r = 10 K; models 3 & 4) 

We now repeat the simulations preformed in models 1 and 
2 in our next two models, 3 and 4 respectively, but with one 
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Fig. 3. A time (measured in Myr) sequence of column density plots showing the commencement of gravitational fragmentation, 
then clumping and finally the fragmented planar gas slab in model 3. 

change. We now hold the post shock gas slab at 10 K. As dis- 
cussed above, we intend to explore the effect of cooling on the 
- dynamical evolution of the gas slab. We saw in models 1 and 2 

above, that the gas slab at 54 K was not susceptible to the gravi- 
_ tational instability. 

However, in model 3, where the MCs collide head-on, we 
observe that the post-collision gas slab undergoes gravitational 
fragmentation. In Fig. 3, we present a time sequence of the gas 
slab after its formation, as seen face-on. As the gravitational in- 
stability starts growing, the slab shows signs of flocculation (f = 
0.65 Myr, t - 0.7 Myr in Fig. 3), and eventually clumps in the 
slab condense out. Some of the smaller clumps merge to form 
larger clumps while some others lie embedded in filamentary 
structures. Most of the clumps being spheroidal, appear non- 
circular in the projection maps (see Fig. 4 below). 

The gravitational instability grows on a time scale much 
smaller than the free fall time of the individual precollision 
clouds, or even the cloud crushing time, t cr Q- We now compare 
observations from this model with corresponding analytic pre- 
dictions. In our simulation we observe that the clump formation 
Fig. 5. A column density plot showing a closer view of the cen- commences after about 0.65 Myr. This is the time required for 
tral region of the fragmented slab in model 3. The (three) black 
dots are the (three) sinks formed in various clumps. 




1 to. = —7^, where v is the precollision speed of a cloud and R c u is 
its radius. 
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the growth of unstable mode, t grawt h, leading to condensation of 
a clump and defined as 

/ L \i 

t ^~[ Go . bfer _ a 2 L -ij • W 

Whitworthefa/(1994). 

Here L, cri aver , and, a/ fl , r r, respectively are the size of the 
clump, average surface density of the slab, and the average sound 
speed in the slab. Plugging in the appropriate values from the 
simulation we get t growt h = 0.6 Myr, which is in agreement with 
the epoch observed in the simulation. We again note that, clump 
formation commences after about 0.65 Myr. Thereafter the grav- 
itational instability grows quickly and the slab fragments into a 
number of clumps. The minimum mass of a fragment is given by 

tl w , .1 layer 

Mfrag = M 2 — — r , (5) 

(G i player) 1 

Whitworth et al (1994). Here pi ayer is the average density of the 
layer and M is the precollision Mach number. Using equation 
(5), we calculate the mass of this fragment, Mf rag , as 0.32 M 
(0.22 M Q ) for the present simulation. The number in brackets 
is the minimum mass of the fragment that condenses out of the 
slab. 

The fragmentation of this slab results in the formation of 36 
big and small clumps. Fig. 4 shows a close-up of the fragmented 
slab. At this epoch a few clumps have already become self gravi- 
tating while a few others, oscillate for some time before collaps- 
ing under self gravity. The column density plot in Fig. 5 shows 



a closer view of the central region of this slab. Three sinks have 
been marked by black dots in as many clumps. Yet some other 
small clumps merge with similar clumps and form larger clumps 
and filaments. To facilitate identification of structure in the slab, 
we have overlayed density contours on the column density plot 
in this figure. 

The slab as seen in Figs. 4 and 5 is fairly evolved and at this 
epoch, the gravitational instability has fully grown. Almost all 
of the gas in the slab ends up in density structures. The aver- 
age density of the structure in the slab is ~10~ 18 g cirT 3 , which 
then suggests that the growth period of the instability as defined 
by Elmegreen (1989), and quoted above, is ~ 0.8 Myr. This is 
comparable to the epoch (0.85 Myr) by which the instability has 
fully grown. 

According to Nakamura, Hanawa & Nakano (1993), the 
characteristic length scale of a filamentary structure is 

H = 0.0035f a ' ayer , Y h 'T T 'pc, (6) 
\0.3kms- 1 A2xl0 6 cm-V V 

where all the symbols have their usual meanings. The average 
temperature of the densest regions in the slab is about 10K while 
"■layer ~ 10 7 cirr 3 . Substituting these values in equation (6) 
above, we get H ~ 0.13 pc. The gravitational fragmentation of 
the filaments produce clumps and the separation between indi- 
vidual clumps is of the order of Jeans length, A j, which for a 
filament is Aj ~ 22H (Nakamura, Hanawa & Nakano 1993). 
In the present case, Aj ~ 0.3 pc. This agrees with the average 
separation between clumps embedded in filamentary structures, 
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Fig. 6. A time (measured in Myr) sequence of column density plots of the oblique slab in model 4. A face-on view of the slab has 
been shown in these plots. We can see that as the slab grows old, it shows signs of fiocculation (t = 0.82 Myr) and subsequently 
fragments (right hand plot in the third row, t — 1 .6 Myrs). The last plot shows a network of filaments and clumps embedded in them. 



as seen in Fig. (3) above. This simulation was terminated after 
t ~ 0.87 Myr by which time seven sinks had formed in various 
clumps. 

Next, in model 4 the cloud collision produces an oblique 
gas slab confined by ram pressure. Like in the previous case 



(model 3), this slab also undergoes gravitational fragmentation. 
However, the fragmentation observed here is somewhat different 
from that reported by Chapman et al (1992) and Pongracic et 
al (1992). According to these authors, the oblique gas slab tum- 
bles about a direction perpendicular to the plane of collision and 
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breaks in to two blobs, which then undergoes secondary frag- 
mentation to form multiples. Such a fragmentation resembles the 
bar mode instability. In the present case, although the oblique 
slab exhibits a similar tumbling motion the final outcome is dif- 
ferent from that reported by the respective authors. 

A possible reason for this observed difference could be, a 
smaller Jeans length in the present case. Our pre-collision clouds 
being more massive produce a slab with greater column density, 
which shortens the Jeans length and enhances its growth rate. In 
this eventuality, the gravitational instability will dominate over 
other unstable modes of the tumbling slab, notably the bar mode. 
The fragmentation of the slab after its formation can be seen in 
Fig. 6. This is a time sequence of column density plots of the gas 
slab, as seen face-on. The nature of fragmentation observed here 
is in consonance with the predictions of Whitworth et al (1994) 
for an off-centre cloud collision, with a small impact parameter. 

As in the previous case, we calculated the dynamical proper- 
ties for this slab, using equations (3), (4) and (5). The length of 
the fastest growing unstable mode in the slab, Lf ast , defined by 
equation (3) above is, indeed much smaller than the radial extent 
of the slab. Clumping in the slab commences after ~ 0.4 Myr, 
while the smallest clump has mass ~ 12 M Q and size ~ 0.1 pc. 
Corresponding analytic calculations yield 0.4 Myr, 1 1 M Q and 
0.12 pc, respectively. The calculated values thus, seem to be in 
good agreement with those observed in the simulation. 

Once the slab becomes unstable, it develops floccules and 
soon fragments. The final state of the slab, post-fragmentation, 
can be seen in Fig. 7 which is a column density plot of the slab 
as seen face-on. The clumps and filaments in it have density of 
order 10~ 19 g cm 4 . Density contours overlayed on this plot elu- 
cidate the structure in it. The average smoothing length of SPH 
particles in these simulations is at least an order of magnitude 
smaller than Lf ast . We therefore think that the observed frag- 
mentation of gas slabs is likely to be physical. 

We find that fragmentation of a cold gas slab leads to clump 
formation and self gravitating clumps may spawn star clusters, 
a view that as also been advanced by numerous workers be- 
fore. See for instance Chapman et al (1992), Bhattal et al (1994) 
and Clarke (1999), among others. Here we can see that frag- 
menting gas slabs also produce dense filaments. In fact Burkert 
& Hartmann (2004) have performed 2-d grid simulations of 
self gravitating, finite gas sheets having different geometries. It 
makes an interesting comparison between two of their models 
viz. the 'static circular sheet' and the 'rotating elliptical sheet', 
with the respective gas slabs in models 3 and 4 discussed in this 
section. 

While the formation of the outer ring like structure in the 
gas slab of model 3 (see Fig. 4), looks similar to the structure 
in their collapsing circular sheet, but unlike the global collapse 
observed in their model the gas slab discussed here undergoes 
gravitational fragmentation. However, the reason behind forma- 
tion of this ring, that of edge material initially piling up, might 
still be valid. Similarly, the Burkert & Hartmann (2004) model 
of collapsing elliptical gas sheet leading to the formation of a 
filament is at variance with the result of model 4, discussed here. 
This is obviously due to the faster growth rate of the gravitational 
instability, as compared to the global collapse time of the gas 
slab. However, the shocked slab (both planar and oblique), where 
the gravitational instability is apparently suppressed by hydro- 
dynamical instabilities, undergoes global collapse (see paper I). 
Hartmann & Burkert (2007) have proposed a rotating ellipti- 
cal sheet model to explain the formation of the Orion Integral 
Filament. 



Observational evidence for similar occurrences have been 
reported, for instance by Lada, Alves & Lada (1999) and 
Elmegreen (2002) about a star cluster embedded in an elongated 
region, IC546. This cluster is thought to have been formed due 
to some dynamical trigger 

Slow cloud collision 

Finally, we examine the case of a very low velocity head- 
on cloud collision. As discussed above, an adiabatic EOS, given 
by equation (2) with y = | is used to model the weak shock 
in this case. Fig. 8 shows a time sequence of column density 
plots of the colliding clouds. Post-collision, the clouds merge 
and form a prolate gas cloud (t = 1 .42 Myrs in Fig. 8). The cloud 
is thermally supported against self gravity. The thermal pressure 
gradient within the cloud is non-isotropic due to its spheroidal 
shape. 

Post-collision, as material from the colliding clouds streams 
in to the gas cloud, it progressively becomes denser and there- 
fore warmer, as dictated by the stiff EOS employed in this model. 
This builds up a thermal reservoir within it, and the cloud, there- 
fore starts re-expanding along the collision axis (f = 1 .96 Myrs, 
t = 2.62 Myr and t = 2.96 Myrs in Fig. 8). During this process it 
flattens, and the re-expansion continues till the thermal pressure 
along the y-axis, arrests the lateral collapse. This observation in 
this simulation is in consonance with the qualitative predictions 
of Ramsay (1961) and Mestel (1965). Note that, the cloud is not 
self gravitating and the lateral collapse is just a phenomenon in- 
volving gravo-thermal balance. 

From the snapshots in Fig. 8, we can also see that some ma- 
terial is ejected, also called jets, from the top and bottom ends of 
the cloud. The qualitative feature of this model is similar to that 
of models 1 and 2 discussed above. However, the post collision 
gas cloud evolves on a much larger timescale than in cases 1 and 
2. This is yet another scenario that could lead to the formation 
of diffuse clouds in the ISM. Star formation can commence in 
such clouds only when they become sufficiently dense to support 
gravitational instability. 

5. Conclusions 

Colliding clouds dissipate kinetic energy and produce gas slabs 
which may fragment, leading to formation of density structures 
of various shapes and sizes. Some of these density concentra- 
tions may collapse and form stars while others may suffer tidal 
disruption. This model does not require injection of turbulence 
to trigger star formation and, stellar feedback from one episode 
of star formation may lead to other such events in the neighbour- 
hood. 

Numerical simulations by Clarke & Gittins (2006) for in- 
stance, have shown that a burst event creates local perturba- 
tions in the galactic disk which in turn generate spiral patterns 
in it. Interference of many such perturbations may create com- 
plex structure, that might evolve on the timescale of a rotation 
period. The evolution of a galaxy essentially depends on the 
global star formation rate. Star burst events also play an impor- 
tant role in the chemistry of the intercloud medium. Powerful 
winds from young star clusters drive shock fronts and, while 
the abundance of molecular species may suffer attrition follow- 
ing the post shock ionisation, the pool of heavier elements will 
be replenished via stellar evolution. Cloud collision is a promi- 
nent mechanism of triggered star formation and supposedly, also 
plays a pivotal role in the distribution of prestellar core masses 
and stellar population in the ISM (McKee 1999; Tan 2000). 

In the present work, having explored the paradigm of low ve- 
locity cloud collisions, we have also tried to ascertain the impor- 
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Fig. 7. A column density plot (f = 1.4 Myrs) showing a close view of the oblique slab in model 4. Fragmentation of the slab leads 
to formation of number of clumps and elongated structures as can be seen with the aid of density contours overlayed on the column 
density plot. The circled region shows a fragmenting clump. 



tance post-shock cooling albeit in a rather crude manner, upon 
the dynamical evolution of the pressure confined slab. Our sim- 
ple investigation shows that, under suitable conditions, gas slabs 
become Jeans unstable and fragment to produce clumps and fila- 
ments. We conclude that colder slabs are more likely to fragment 
otherwise, the slab simply re-expands and ends up as a diffuse 
cloud in dynamical equilibrium. The former scenario is more in- 
teresting from the perspective of star formation, either multiples 
or larger A^-body clusters. Density perturbations in the gas slab 
grow purely from white noise and therefore do not require any 
external trigger. In the latter scenario however, the diffuse gas 
clouds may become self gravitating, when crushed by sufficient 
external pressure else, clouds of various shapes simply lie scat- 
tered in the ISM. 

Models 1, 2, and 5 in the present work, belong to this lat- 
ter paradigm while models 3 and 4, to the former. Figs. 4 and 7 
respectively, show the fragmented gas slabs in models 3 and 4. 
Clumps and well defined, elongated density structures are evi- 
dent in both. Indeed, a few clumps in either slabs, have become 
self gravitating by the time respective calculations were termi- 
nated. The collapse of individual clumps however, could not be 
followed further due to shortage of resources. However, frag- 
mentation of gas slabs is evidently a propitious mode of clump 
formation. Clumps may eventually spawn stars and larger clouds 
may become wombs for larger star clusters. 

The results of our off-centre cloud collision experiments 
(both models 2 and 4) differ significantly from similar exper- 
iments performed by Lattanzio et al (1985), Chapman et al 



(1992) and Bhattal et al (1998) among others. These authors 
report a bar mode fragmentation of the post-collision oblique 
pressure compressed slab, which is different from the fragmenta- 
tion observed in our model. On the contrary, we observe that the 
oblique gas slab either, simply re-expands to finally produce a 
diffuse gas cloud (model 2) or, fragments gravitationally to form 
clumps and filaments (model 4). We attribute this difference to 
the fact that the mass of the slab which, in our simulation is at 
least an order of magnitude greater than that in the work cited 
above. As a result of larger mass and therefore higher column 
density, the length of the unstable mode is shortened and grows 
much faster. 

We admit that the investigation presented here is too simpli- 
fied and dwells purely on a gravo-thermal treatment of the prob- 
lem. At best, we have a skewed picture before us for a more elab- 
orate future study. For an elaborate study of the subject, mag- 
netic field needs to be included and post-shock radiative cooling 
needs to be better treated (c. / Vazquez-Semadeni et al 2007). 
Also, collapsing clumps can be followed further by employing 
a radiative transfer scheme (c. / Stamatellos et al 2007), rather 
than a simple barytropic EOS used here. 
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Fig. 8. A time (measured in Myrs) sequence of column density plots in model 5. Collision of clouds, followed by the formation of 
a prolate gas cloud (f = 1.07 Myrs) and then its re-expansion is evident from these plots. The expansion of the post-collision gas 
cloud terminates in a diffuse cloud (f = 2.96 Myrs). Material squirting from the top and bottom ends (jets) of the slab are also visible 
in these plots. 
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